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Abstract. We present a multiwavelength study of the massive star forming region associated with IRAS 
06055+2039 which reveals an interesting scenario of this complex where regions are at different stages of evolution 
of star formation. Narrow band near-infrared (NIR) observations were carried out with UKIRT-UFTI in molec- 
ular hydrogen and Br7 lines to trace the shocked and ionized gases respectively. We have used 2MASS JHKs 
data to study the nature of the embedded cluster associated with IRAS 06055-1-2039. We obtain a power-law 
slope of 0.43±0.09 for the i^s-band Luminosity Function (KLF) which is in good agreement with other young 
embedded clusters. We estimate an age of 2 - 3 Myr for this cluster. The radio emission from the ionized gas 
has been mapped at 610 and 1280 MHz using the Giant Metrewave Radio Telescope (GMRT), India. Apart from 
the diffuse emission, the high resolution 1280 MHz map also shows the presence of several discrete sources which 
possibly represent high density clumps. The morphology of shocked molecular hydrogen forms an arc towards 
the N-E of the central IRAS point source and envelopes the radio emission. Submillimetre emission using JCMT- 
SCUBA show the presence of a dense cloud core which is probably at an earlier evolutionary stage compared 
to the ionized region with shocked molecular gas lying in between the two. Emission from warm dust and the 
Unidentified Infrared Bands (UIBs) have been estimated using the mid-infrared (8 - 21 /im) data from the MSX 
survey. From the submillimetre emission at 450 and 850 nm the total mass of the cloud is estimated to be ~ 7000 
- 9000 M0. 

Key words, infrared: ISM - radio continuum: ISM - ISM: H II regions - ISM: individual objects: IRAS 06055-1-2039 



1. Introduction 

Massive stars are preferentially formed in dense cores of 
molecular clouds. They remain deeply embedded in the 
prenatal molecular gas and obscuring dust and their pre- 
main sequence time scales are much shorter compared to 
the low mass stars. The luminous high mass stars also 
affect the parent cloud. In addition, massive stars do 
not form in isolation but often in clusters and associa- 
tions. All these factors contribute in making the study of 
the formation mechanisms of these systems very difficult. 
Multiwavelength studies, therefore, hold the potential to 
probe these complexes at different depths and unravel the 
least understood aspects of massive star formation. IRAS 
06055-1-2039 (G189.78-f0.34, RAFGL 5179) is a massive 
star forming region chosen from the catalog of massive 
young stellar objects by Chan ct al. (1996). G189. 78+0.34 
is listed as an ultracompact (UC) HII region (Shepherd & 
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Churchwell 1996; Bronfman et al. 1996). It belongs to the 
Gem OBI molecular cloud complex and is a part of the 
extended HII region Sh 252. It is associated with S252 A 
which is one of the six compact radio sources in Sh 252 
revealed from the 5 GHz aperture synthesis observations 
by FeUi et al. (1977). IRAS fluxes yield a far infrared lu- 
minosity of ^ lO^'Lo (Carpenter et al. 1995) for this star 
forming region. There are several kinematic distance esti- 
mates used in literature for this source which range from 
1.25 kpc (Mirabel et al. 1987) to 2.9 kpc (Wouterloot & 
Brand 1989). In this paper we use the value of 2.6 kpc 
(Wu et al. 2001) which is the most widely used distance 
estimate for this source. 

This high mass star forming region has been observed 
as part of many surveys. H2O maser (Kompe et al. 1989; 
Lada & Wooden 1979) and the 6.7 GHz methanol maser 
(Szymczak et al. 2000a) have been detected towards IRAS 
06055+2039. Positive detection has also been made in 
SiO (Harju et al. 1998), CO (2 - 1) (Wu et al. 2001) 
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and CS (2 - 1) (Bronfman et al. 1996). Zinchenko et al. 
(1998) in their study of dense molecular cores have also 
mapped this source in CS (2 - 1). CO maps of Shepherd 
& Churchwell (1996) do not show evidence of any bipo- 
lar outflows. Search for the 6 cm (Szymczak et al. 2000b) 
and 5 cm (Baudry et al. 1997) OH masers show nega- 
tive results. The IRAS low resolution spectra show a rela- 
tively red continuum from 13 to 23fim with the presence 
of an emission feature at 11.3 /xm which is attributed to 
the presence of Polycyclic Aromatic Hydrocarbon (PAH) 
molecules (Kwok et al. 1997). 

In this paper, we present a multiwavelength study of 
this star forming region. In Sect. [2 we present the narrow- 
band near-infrared (NIR) and radio continuum observa- 
tions and a brief description of the related data reduction 
procedures. In Sect.|Hl we discuss other available datasets 
used in the present study. Section 0] gives a comprehen- 
sive discussion on the results obtained and in Sect.[5l we 
summarize the results. 



2. Observations and Data Reduction 

2.1. Near-Infrared Observations 

Narrow-band NIR observations were carried out in the 
rotational- vibrational line of molecular hydrogen - H2 (1- 
0)S1 (A=2.12/im, AA=0.02/im), hydrogen recombination 
line of Br7 (A=2.16/im, AA=0.02/im) and K continuum 
(A=2.27/im, AA=0.02;um) at the 3.8m United Kingdom 
Infrared Telescope (UKIRT^), Hawaii. Observations were 
carried out on 16 Feb 2002 under the UKIRT Service 
Programme (Proposal No. 1459). The instrument used 
was the UKIRT Fast Track Imager (UFTI), which is a 
1 - 2.5 nm camera with a 1024x1024 HgCdTe array hav- 
ing a plate scale of 0.091" per pixel. The 'JITTER-SELF- 
FLAT' Data Reduction (DR) recipe was used. This script 
takes imaging observations comprising of 9 jittered object 
frames and a dark frame. Flat field is then created from 
the sequence of jittered object frames. The final image is a 
mosaic generated from the 9 frames after dark subtraction 
and flat flelding and has a total field of view of ~ 4x4 
arcmin^. For our observations, the integration time was 
100s in each band. To obtain pure emission line images 
it is essential to subtract out the contribution from the 
continuum. This is done by subtracting the K continuum 
image from the H2 and the Br7 images after proper align- 
ment and PSF matching. Figure ^ shows the continuum 
subtracted H2 and Br7 images of the central field around 
IRAS 06055+2039. The continuum subtracted H2 image 
displays a prominent arc towards the N-E of the central 
source. The continuum subtracted Br7 image shows the 
presence of faint diffuse emission surrounding the central 
bright source. 



^ The United Kingdom Infrared Telescope is operated by the 
Joint Astronomy Centre on behalf of the U.K. Particle Physics 
and Astronomy Research Council. 



Table 1. Details of the radio interferometric continuum 
observations. 



Details 


1280 MHz 


610 MHz 


Date of Obs. 


3 Aug 2003 


21 May 2002 


Primary beam 


26' 


54' 


Cont. bandwidth (MHz) 


16 


16 


Synth, beam 


4".0 X 4".0 


6". 3x4". 3 


Position angle, (deg) 




58 


Peak Flux (mjy/beam) 


5.7 


7.9 


rms noise (mJy/beam) 


0.3 


0.4 



2.2. Radio Continuum Observations 

In order to probe the ionized gas component, radio contin- 
uum interferometric mapping of the region around IRAS 
06055+2039 was carried out using the Giant Metrewave 
Radio Telescope (GMRT) array, India. The GMRT has a 
"Y" shaped hybrid configuration of 30 antennas, each of 
45 m diameter. There are six antennas along each of the 
three arms (with arm length of ~ 14 km). These provide 
high angular resolution (longest baseline ~ 25 km). The 
rest of the twelve antennas are located in a random and 
compact 1x1 km^ arrangement near the centre and is sen- 
sitive to large scale diffuse emission (shortest baseline ^ 
100 m). Details of the GMRT antennae and their configu- 
rations can be found in Swarup et al. (1991). Observations 
were carried out at 1280 and 610 MHz. The radio sources 
3C48 and 3C147 were used as the primary flux calibrators 
while 0532+194 and 0432+416 were used as phase calibra- 
tors for the 1280 and 610 MHz observations respectively. 

Data reduction was done using AIPS. The data sets 
were carefully checked using tasks UVPLT and VPLOT 
for bad data (owing to dead antennae, bad baselines, in- 
terference, spikes etc) . Subsequent editing was carried out 
using the tasks UVFLG and TVFLG. Maps of the field 
were generated by Fourier inversion and cleaning using 
the task IMAGR. Several iterations of self calibration were 
carried out to obtain improved maps. 

Figure 121 shows the radio continuum images at 1280 
and 610 MHz with synthesized beam sizes of 4".0 x 4".0 
and 6". 3x4". 3 respectively. The rms noises in the maps are 
0.3 (1280 MHz) and 0.4 (610 MHz) mJy/beam. Table[l] 
gives the details of the observations and the maps. 

3. Otiier Available Datasets 

3.1. Near-Infrared Data from 2MASS 

NIR (JHKs) data for point sources around IRAS 
06055+2039 have been obtained from the Two Micron 
All Sky Survey^ (2MASS) Point Source Catalog (PSC). 
Source selection was based on the 'read-flag' which gives 

^ This publication makes use of data products from the 
Two Micron All Sky Survey, which is a joint project of the 
University of Massachusetts and the Infrared Processing and 
Analysis Center/California Institute of Technology, funded by 
the NASA and the NSF. 
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Fig. 1. The continuum subtracted Ha (l-O)Sl (left) and Br7 (right) images of the field around IRAS 06055+2039. 
The images were smoothed to a resolution of 0.3" x 0.3" (3x3 box) to improve the contrast of the faint features. 
Residuals of continuum subtraction around stars are seen in both the images. The plus sign and the asterisk in both 
the images mark the position of the IRAS point source and the central brightest star in the NIR cluster respectively. 
Artifacts (possibly ghost images of the central bright star) are seen in II2, Br7 and the K continuum frames. The 
resulting residuals in both the continuum subtracted images are indicated with lines. 



the uncertainties in the magnitudes. In our sample we re- 
tain only those sources for which the 'read-flag' values are 
1-6. The 2MASS data have been used to study the em- 
bedded cluster associated with IRAS 06055+2039. 

3.2. Mid-Infrared Data from MSX 

The Midcourse Space Experiment'^ (MSX) surveyed the 
Galactic plane in four mid infrared bands - A (centre 
A: 8.28 /xm; 50% peak intensity range: 6.8 ~ 10.8 /xm), C 
(12.13Aim; 11.1 - 13.2/xm), D (14.65^m; 13.5 - 15.9/im) 
and E (21.34/im; 18.2 - 25.1/im) at a spatial resolution 
of ~ 18" (Price et al. 2001). Two of these bands (A & 
C) cover the Unidentified Infrared Bands (UIBs) at 6.2, 
7.7, 8.7, 11.3 and 12.7/im. The integrated flux densities 
of IRAS 06055+2039 in these bands are listed in Table El 
The MSX images in these four bands for the region around 
IRAS 06055+2039 have been used to study the emission 
from the UIBs and to estimate the spatial distribution of 
temperature and optical depth of the warm interstellar 
dust. 



^ This research made use of data products from the 
Midcourse Space Experiment. Processing of the data was 
funded by the Ballistic Missile Defense Organization with ad- 
ditional support from NASA Office of Space Science. This re- 
search has also made use of the NASA/ IPAC Infrared Science 
Archive, which is operated by the Jet Propulsion Laboratory, 
Caltech, under contract with the NASA. 



3.3. Mid- and Far-infrared Data from IRAS 

The data from the IRAS survey in the four bands (12, 25, 
60 and 100 /im) for the region around IRAS 06055+2039 
were HIRES processed (Aumann et al. 1990) to obtain 
high angular resolution maps. These maps have been used 
to determine the spatial distribution of dust colour tem- 
perature and optical depth. The integrated flux densities 
from the HIRES processed images and IRAS-PSC are also 
given in Table El 

3.4. Sub-mm Data from JCMT 

Submillimetre observations using the Submillimetre 
Common-User Bolometer Array (SCUBA) instrument of 
the James Clerk Maxwell Telescope'' (JCMT) were re- 
trieved from the JCMT archives and processed using 
their standard pipeline SCUBA User Reduction Facility 
(SURF). JCMT-SCUBA observations for the data used 
in our study were carried out on 25 Oct 2000. Uranus was 
used as primary flux calibrator for the maps. FigureOdis- 
plays the spatial distribution of cold dust emission at 450 
and 850 /xm. Large atmospheric extinction correction has 
been applied to obtain the 450 /xm map. 



^ This paper makes use of data from the James Clerk 
Maxwell Telescope Archive. The JCMT is operated by the 
Joint Astronomy Centre on behalf of the UK particle Physics 
and Astronomy Research Council, the National Research 
Council of Canada and the Netherlands Organisation for Pure 
Research. 
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32 31 
RIGHT ASCENSION (J2000) 
Cont peak flux = S.7307E-03 JY/BEAM 
Levs = 3.000E-04 " (2.500, 3.500, 5, 7, 9, 10, 11, 
13, 14, 15, 16, 17, 17.50, 18.50) 




32 31 
RIGHT ASCENSION (J2000) 
Cant peak flux = 7.9258E-03 JY/BEAM 
Levs = 4.000E-04 • (2.500, 3.500, 5, 7, 9, 1 1 , 13, 
15, 16, 16.60, 17.50, 18, 18.50, 19.50) 



Fig. 2. High resolution radio continuum maps at 1280 MHz (left) and 610 MHz (right) for the region around IRAS 

06055+2039. The numbers Rl, R2, R7 mark the position of discrete points sources (see text). The rms noises at 

1280 and 610 MHz are ~ 0.3 and 0.4 mJy/beam respectively. The plus sign in both the images marks the IRAS point 
source position. 

Table 2. Infrared flux densities for IRAS 06055+2039 





MSX" 


images 






IRAS 


imaj 


;es 




Wavelength (/im) 


8.3 12.1 


14.7 


21.3 




12 


25 


60 


100 


Flux density (Jy) 


26 34 


24 


68 


PSC 


15.6 


76.8 


1032 


1715 










HIRES" 


38 


113 


1280 


1691 



" Flux densities were obtained by integrating over a circular region of diameter 3' centered on the peak. 



4. Results and Discussion 

4.1. Embedded Cluster in the Near-Infrared 

4.1.1. Radial Profile and Stellar Surface Number 
Density 

The 2MASS Kg - band image of the region around IRAS 
06055+2039 is shown in Figure^ We see the presence of 
a diffuse emission region harbouring an infrared cluster. 

This cluster is also listed in the catalog of embed- 
ded infrared clusters compiled by Bica et al. (2003). We 
use the 2MASS JHKg data to study the nature of this 
cluster. To determine the Kg - band radial profile and 
the stellar surface number density (SSND), we have se- 
lected sources which are detected in the Kg band. To 
estimate the cluster radius we select a large region of 
radius 300" centered on IRAS 06055+2039 (02000.0 = 
06'^08'"32M; Jzooo.o = +20°39'18"). To account for the 
contribution from the field stars we select a control field 
(a2ooo.o = 06h09™52^0 ; Jzooo.o = +20°39'18") which is ~ 
20' to the east of IRAS 06055+2039. FigureEl shows the 



radial profile of the stellar density in log-log scale. This 
profile was created by counting the number of stars in 10" 
annuli and normalizing by the annulus area. We fit two 
models to the surface density radial profile - the King's 
profile and the inverse radius {r~^) model. Neglecting the 
tidal radius, the King's profile can be written as 



/(r) = a + 



fo 



1 + (r/re)2 



(1) 



where, /o is the core concentration at radius zero, rc is 
the core radius and a is a constant for the background 
offset. As seen in Fig. El both the models describe the 
density distribution fairly well. However, the model 
has a better overall fit (reduced = 1-2) as compared to 
the King's model (reduced = 3.3). Several studies have 
shown that young embedded clusters have been fitted by 
r^^ profiles (e.g McCaughrean & Stauffer 1994; Lada & 
Lada 1995) or, by both the inverse radius model and the 
King's model (e.g Horner et al. 1997; Teixeira et al. 2004; 
Baba et al. 2004). As discussed in Baba et al. (2004), the 
dependence is likely to be a reminiscent of the parental 
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Fig. 3. Contour maps showing the spatial distribution of dust emission at 450 /im (left) and 850 /im (right) for the 
region around IRAS 06055+2039. The contour levels are at 5, 7, 9, 10, 15, 20, 30, 40, 50, 60, 70, 80 and 90% of the 
peak value of 29.45 and 6.83 Jy/beam at 450 and 850 /xm respectively. The FWHMs of the symmetric 2-D Gaussian 
beams are 7". 8 and 15". 2 for the two wave bands. The plus sign in both the images marks the position of the IRAS 
point source. 
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Fig. 4. The 2MASS Kg - band image of the region around 
IRAS 06055+2039. The presence of a diffuse emission re- 
gion and a NIR cluster is seen. The plus sign marks the 
position of the IRAS point source. 



cloud core whereas the King's profile represents systems in 
dynamical equilibrium. Hence, the better fitting of the 
model could possibly suggest that the cluster associated 
with IRAS 06055+2039 is not yet in complete dynamical 
equilibrium. 



Within errors, the cluster profile merges with the field 
star level at 85" which translates to ~ 1.1 pc at a dis- 
tance of 2.6 kpc. We take this as the cluster radius. The 
background level as estimated from the control field is 
^ 9 stars pc^^ which is in reasonable agreement with the 
value of 11.8 ± 1.6 stars pc^^ yielded by the King's profile 
fitting. The King's profile fitting also gives a core radius 
rc ^ 0.1 pc. The core radius is a scale parameter and 
depends mainly on cluster parameters like density, lumi- 
nosity, total mass etc. Several stiidics have shown that 
the core radius of the cluster is also correlated with its 
age. In their study of young clusters in the LMC, Elson 
et al. (1989) show that the core radii increases between 
^ 10^ and lO'' yr and then begin to decrease again. Such 
trends in core radius evolution has also been discussed by 
Wilkinson et al. (2003) for the LMC clusters and Mackey 
& Gilmore (2003) for clusters in the SMC. These authors 
have shown that apart from the general increasing trend, 
the spread in the core radii also increases with the age of 
the cluster. Teixeira et al. (2004) and Baba et al. (2004) 
derive core radius values of 0.05 and 0.08 pc for the clus- 
ters NGC 2316 and RCW 36 respectively The age esti- 
mate for both these clusters is 2 - 3 Myr. For a compara- 
tively older (5 - 10 Myr) cluster NGC 2282, Horner et al. 
(1997) obtain a core radius of 0.19 pc. The above values 
of core radii seem to suggest that the cluster associated 
with IRAS 06055+2039 is likely to be of the same age (~ 
2-3 Myr) as NGC 2316 and RCW 36. 

The number of stars detected in the Kg band 

within 85" radius of the cluster is 114. The total 
background population is 34 [= 9(background) x tt x 
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Fig. 5. The radial profile of the surface number density 
for the cluster associated with IRAS 06055+2039 in log- 
log scale. Also plotted are the two fitted models - King's 
model (solid) and inverse radius (dashed). As compared 
to the King's model, the inverse radius model fits the ra- 
dial profile much better. The horizontal dotted line cor- 
responds to the background field star level which is 9 
stars pc"'^. Statistical errors are shown. 

(l.l)^(areaof cluster)]. Hence, the total number of cluster 
members is estimated to be 80. This yields a space den- 
sity of ~ 14 stars pc~'^. However, it should be noted here 
that this density is a lower limit as we are not completely 
sampling the stellar population at fainter magnitudes due 
to the 2MASS sensitivity limits. 

Figure El shows the surface density map of the region 
around IRAS 06055+2039. We see the presence of two 
peaks. The prominent peak which lies to the west coin- 
cides spatially with the central bright source. The sec- 
ondary peak, which lies ^ 50" to the east of the main 
peak, coincides well with the departure seen in the sur- 
face density profile from the model profiles around a ra- 
dial distance of ~ 0.6 pc. This secondary peak is situated 
close to the edge of the sub-mm peaks presented in Fig. O 
The stellar surface density distribution exhibits a centrally 
condensed-type structure rather than a hierarchical-type 
structure (Lada & Lada 2003). This is consistent with the 
fact that the derived radial profile fits reasonably well both 
to the King's model as well as the inverse radius model. 
Hierarchical-type complexes do not follow any well de- 
fined profile. Though at a less significant level, centrally 
condensed clusters have also been seen to show the pres- 
ence of structures. For example, Lada & Lada (1995) show 
the presence of satellite subclusters in the outer regions of 
IC 348. This is consistent with the double peaked struc- 
ture that we see for the cluster associated with IRAS 
06055+2039. 

4.1.2. Colour-Colour (CC) and Colour-Magnitude 
(CM) Diagrams 

The {H — K) versus (J — H) CC diagram for the cluster 
associated with IRAS 06055+2039 is shown in Fig. In 
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Fig. 6. The contour map of the stellar surface number 
density obtained by counting stars in a 10" x 10" ( ^ 
0.1 X 0.1 pc) grid for the cluster associated with IRAS 
06055+2039. The contours are from 30 to 140 stars pc"^ 
in steps of 20 stars pc~^. The lowest contour is at the ia 
level. The positions of the IRAS point source (plus), the 
radio peak (open triangle) and the sub-mm peak (star) 
are shown in the figure. 

this figure we have plotted 53 sources which have good 
quality photometric magnitudes in all the three JHKg 
bands (2MASS 'read-flag' value of 1 - 3). Henceforth, for 
the analysis of 2MASS data, we will be using only good 
quality photometric data which have the above 'read-flag' 
values. For clarity we have classified the CC diagram into 
three regions (e.g. Sugitani et al. 2002; Ojha et al. 2004a 
& b). The "F" sources are located within the reddening 
bands of the main sequence and the giant stars. These 
sources are generally considered to be either field stars. 
Class HI objects or Class II objects with small NIR excess. 
"T" sources populate the region redward of the "F" region 
but blueward of the reddening line corresponding to the 
red end of the T Tauri locus. These sources are classical 
T Tauri stars (Class II objects) with large NIR excess or 
Herbig AeBe stars with small NIR excess. Redward of the 
"T" region is the "P" region which are mostly protostar- 
like Class I objects and Herbig AeBe stars. Majority of 
sources in our sample are almost equally distributed in 
the "F" and the "T" regions, whereas, only four lie in the 
"P" region. A total of 18 out of 80 (22%) sources show in- 
frared excess (i.e. sources populating the "T" and the "P" 
regions). However, it is important to note here that this 
NIR excess fraction is just the lower limit as several clus- 
ter members detected in the Kg band were not detected in 
the other two shorter wavelength bands and hence are not 
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part of this sample. The NIR excess in pre-main sequence 
(PMS) stars is due to the optically thick circumstellar 
disks /envelopes. These disks/envelopes become optically 
thin with age hence the fraction of NIR excess stars de- 
creases with age. For very young (< 1 Myr) embedded 
clusters the fraction is ~ 50 % (Lada et al. 2000; Haisch 
et al. 2000) and it decreases to ~ 20 % for more evolved (2 

- 3 Myr) clusters (Haisch et al. 2001; Texeira et al. 2004). 
The fraction of stars with NIR excess seen in this cluster 
suggests an upper limit 2-3 Myr on the age which is 
reasonably consistent with that suggested from the core 
radius value. The age estimates for the cluster NGC 2175 
associated with the extended HII region Sh 252 are 2 Myr 
(Grasdalen & Carrasco 1975) and 1-2 Myr (Kompe et 
al. 1989) which agrees rather well with our estimates for 
the cluster associated with IRAS 06055+2039 which also 
belongs to the Sh 252 complex. 

We have calculated the extinction by de-reddening the 
stars in the CC diagram. The stars are shifted to a line 
drawn tangential to the turn-off point of the main se- 
quence locus (see Fig.[7||. The amount of shift gives an 
estimate of the extinction of individual stars. The extinc- 
tion values range from Ay ^ to 13 mag with an average 
foreground extinction of Ay ^ 7 mag. The range of ex- 
tinction values obtained shows up as the spread of stars 
along the reddening band in the CC diagram. This indi- 
cates that the cluster is partially embedded (Teixeira et 
al. 2004). Also according to Lada & Lada (2003), low ex- 
tinction values {Ay ~ 1 - 5 mag) are typical of partially 
embedded clusters. 

Figure IHlshows the {H—K) versus K colour-magnitude 
(CM) diagram for 79 sources with good quality HK mag- 
nitudes. Using the zero age main sequence (ZAMS) loci 
and the reddening vectors, we estimate the spectral type 
of the brightest star in the cluster to be ~ BO. 5. This is 
the central IRAS point source and from the NIR estimates 
seems to be the most massive star in the cluster. A simi- 
lar estimate for the spectral type is also obtained from the 
analysis of the (J — H) versus J CM diagram. 

4.1.3. Ks - Band Luminosity Function 

We use the 2MASS Kg - band star counts to derive the Ks 

- band luminosity function (KLF) for the embedded clus- 
ter. In order to obtain the KLF of the cluster it is essential 
to account for the background and foreground source con- 
tamination. For this purpose we use both the Besangon 
Galactic model of stellar population synthesis (Robin et 
al. 2003) and the observed control field star counts. We 
use the same control field as described in Sect 14.1.11 

Star counts were predicted using the Besangon model 
in the direction of the control field. We have checked the 
validity of the simulated model by comparing the model 
KLF with that of the control field and found both the 
KLFs to match rather well. As mentioned in the previous 
section, the average foreground extinction is determined 
to be Ay~ 7 mag. Hence, assuming spherically symmetric 
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Fig. 7. Colour-colour diagram of the infrared cluster in 
the IRAS 06055+2039 region. The two solid curves repre- 
sent the loci of the main sequence (thin line) and the giant 
stars (thicker line) derived from Bessell & Brett (1988). 
The long-dashed line is the classical T Tauri locus from 
Meyer et al. (1997). The parallel dotted hues are reddening 
vectors with the crosses placed along these lines at inter- 
vals corresponding to five magnitudes of visual extinction. 
We have assumed the interstellar reddening law of Rieke 
& Lebofsky (1985) {Aj/Ay = 0.282; Ah /Ay = 0.175 and 
AK/Ay = 0.112). The short-dashed line represents the 
locus of the Herbig AeBe stars (Lada & Adams 1992). 
The plot is classified into three regions namely "F", "T" 
and "P" (see text for details). The colours and the curves 
shown in the figure are all transformed to the Bessell & 
Brett (1988) system. The solid line shown is drawn tan- 
gential to the turn-off point of the main sequence locus. 
The arrow points to the position corresponding to the cen- 
tral brightest star in the cluster. The mean photometric 
errors are shown in the upper right corner. 



geometry, the background population is then seen through 
a cloud with extinction upto Ay^ 14 mag (7 x 2). Model 
simulations with Ay — mag and d < 2.6 kpc gives the 
foreground contamination. The background population is 
generated with Ay = 14 mag and d > 2.6 kpc. We deter- 
mine the fraction of the contaminating stars (foreground 
+ background) over the total model counts. This fraction 
is used to scale the observed control field and subsequently 
the star counts of the modified control field are subtracted 
from the KLF of the cluster to obtain the final corrected 
KLF which is shown in the left panel of Fig. |51 

The right panel of Fig. |2| shows the KLF plotted as 
logN versus the K magnitude. The KLF can be written 
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Fig. 9. Left: The corrected Kg - band luminosity function (KLF) for the cluster around IRAS 06055+2039 is shown 
(solid line). The dotted line is the luminosity function without foreground/background correction. Right: The KLF 
shown as logN versus the Kg magnitude. The straight line is the least squares fit to the data points in the magnitude 
range 11.5 - 14.5. 



as a power-law 



dN{K,) 



cx 10 



(2) 



where, the left side of the equation denotes the number 
of stars per unit magnitude bin and a is the slope of 
the power-law. A linear least squares fitting algorithm is 
used to fit the above power-law to the KLF in the mag- 
nitude range Ks = 11.5 to 14.5. We obtain a value of a 
= 0.43±0.09 for the cluster. The least square fitting is 
done taking into account the statistical errors on individ- 
ual data points. 

Within the quoted errors the estimate of the power-law 
slope is consistent with the average value of slopes (a ~ 
0.4) obtained for other young clusters (Lada et al. 1991; 
Lada & Lada 1995; Lada & Lada 2003). The power-law 
slope values obtained for other embedded clusters like W3 
main and NGC 7538 are significantly lower (a ~ 0.17 to 
0.33 - Ojha et al. 2004a & b; Balog et al. 2004). However, it 
should be noted here that these clusters are much younger 

1 Myr) and the surveys are deeper {Kg < 17.5) which 
probe the low mass stellar population down to ^ O.IMq. 

We estimate the masses of the sources in the cluster by 
comparing them with the evolutionary models of Palla & 
Stabler (1999). Figure Uni shows the (J - H) versus J CM 
diagram for the cluster field. We use the J-band magni- 
tudes rather than Ks because it is less affected by emission 
from circumstellar material. The solid curve represents the 
ZAMS isochrone for a 2 Myr cluster from Palla & Stabler 
(1999). Majority of sources have a typical mass of ~ 2M0 



which we assume to be representative of the stellar pop- 
ulation in the cluster. For a log-normal IMF, the power 
law slope (7), which is 1.35 for a Salpeter IMF, is variable 
and is given by 7 = 0.94 + 0.941og(m^), where m^, is the 
stellar mass (Miller & Scalo 1979; Lada et al. 1993). Using 
this relation, we estimate 7 ~ 1.2 for this mean mass. The 
power-law slope for the mass to luminosity relation is f3 
« 1 for clusters of age ^ 10^ Myr (Simon et al. 1992; 
Lada et al. 1993). The corresponding slope of the KLF, 
a (= 7/2.5/3), is 0.48 which is consistent with the value 
obtained from the least squares fit to the KLF. As is seen 
from Fig. 1101 the mass of the majority of sources in the 
cluster are below 2.5 M© and the lowest mass limit is 
0.4 M0 from our sample. 

4.2. Spatial Distribution of UIBs from MSX Data 

We have used the scheme developed by Ghosh & Ojha 
(2002) to extract the contribution of UIBs (due to the 
Poly cyclic Aromatic Hydrocarbons (PAHs)) from the mid- 
infrared images in the four MSX bands. The emission from 
each pixel is assumed to be a combination of two compo- 
nents. The first is the thermal continuum from dust grains 
(gray body) and the second is the emission from the UIB 
features falling within the MSX bands. The scheme as- 
sumes the dust emissivity to follow the power law of the 
form e\ oc and the total radiance due to the UIBs 
in band C to be proportional to that in band A. A self 
consistent non-linear chi-square minimization technique is 
used to estimate the total emission from the UIBs, the 
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Fig. 8. Colour-magnitude diagram of the infrared cluster 
in the IRAS 06055+2039 region. The nearly vertical solid 
lines represent the zero age main sequence (ZAMS) loci 
with 0, 20 and 40 magnitudes of visual extinction cor- 
rected for the distance. The slanting lines show the red- 
dening vectors for each spectral type. The magnitudes and 
the ZAMS loci are all plotted in the Bessell & Brett (1988) 
system. The arrow points to the position corresponding to 
the central brightest star in the cluster. The mean photo- 
metric errors are shown in the upper right corner. 

temperature and the optical depth. The spatial distribu- 
tion of emission in the UIBs with an angular resolution 
of ~ 18" (for the MSX survey) extracted for the region 
around IRAS 06055+2039 is shown in Fig. El 

Comparing the morphology with the radio continuum 
maps, we see that the emission from the UIBs is much 
more extended though the gross morphologies are similar 
with a relatively steep intensity gradient towards the N- 
E and a smoothly decreasing intensity distribution to the 
S-W. The peak position of the UIB distribution matches 
rather well with the radio peak. 

The integrated emission from the region around IRAS 
06055+2039 in the UIB features within the band A of 
MSX (viz., 6.2, 7.7 & 8.6 /xm) is found to be 2.85 xlO^^^ 
Wm~^ (see Fig. Ill|l . For comparison, we have estimated 
the emission in individual UIB features from the IRAS 
LRS spectrum (covering 8 - 22 /zm), for this source. The 
total emission in the 7.7 and 8.6 fim features is ~ 7.33 
xlO-12 Wm-2, in the 11.3^m feature is ~ 1.32 xlO'^^ 
Wm~^ and an upper limit for the 12.7/J,m feature is 5.9 
xlO~^^ Wm~^. The last value is an upper limit due to 
possible contamination from the Ne [II] line at 12.8 /im. 
Hence, the UIB emission extracted from the MSX band 
A is 39% of the estimate from LRS. This is reasonable 
considering the larger effective field of view for the latter. 
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Fig. 10. The mass spectrum for the cluster associated 
with IRAS 06055+2039. The sohd curve represents the 
model isochrone from Palla & Stabler (1999) for a 2 Myr 
PMS stellar population for the mass range 0.1-3 M©. 
The slanted dotted lines are the reddening vectors for 
2.5 and 0.2 Mq PMS stars respectively. Sources which lie 
above the reddening vector for 2.5 Mq are luminous mas- 
sive stars. The arrow points to the position corresponding 
to the central massive star in the cluster. 



4.3. Emission from Ionized Gas 

The radio continuum emission from the ionized gas asso- 
ciated with the region around IRAS 06055+2039 at 1280 
and 610 MHz is shown in Fig. El The integrated flux den- 
sities from these maps are 183 and 282 mJy at 1280 and 
610 MHz respectively. It should be noted here that the 
flux densities are obtained by integrating upto 3cr level, 
where a is the rms noise of the maps. The integrated flux 
densities obtained from the GMRT maps are consistent 
with the results of Felh et al. (1977) and White & Gee 
(1986). FeUi et al. (1977) get values of 230 and 205 mJy 
at 1415 and 4995 MHz respectively. White & Gee (1986) 
estimate the total integrated flux density at 5 GHz to be 
198 mJy. The contour maps display a cometary morphol- 
ogy with a bright arc-shaped edge on the N-E side and a 
smoothly decreasing intensity distribution on the opposite 
side. The cometary morphology is more clearly seen in the 
610 MHz map. Such a morphology implies that the HII 
region is ionization bounded towards the N-E and density 
bounded towards the S-W. The position of the brightest 
radio peak matches with the central bright star seen in 
the infrared images. Similar morphology is seen from the 
4995 MHz map of Felh et al. (1977). They report the pres- 
ence of several dense clumps. According to them the radio 
continuum emission from this HII region (S252 A) is spa- 
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Fig. 11. Spatial distribution of the total radiance in the 
UIBs for the region around IRAS 06055+2039 as extracted 
from the MSX four band images. The contour levels are 
at 5, 10, 20, 30, 40, 50, 60, 70, 80, 90 and 95 % of the peak 
emission of 4.1 x 10~^ Wm~^ Sr~^. The peak position of 
the UIB distribution (filled triangle), the 1280 MHz radio 
emission (open star) and the IRAS point source position 
(plus) are also indicated. 

tially separate from the larger extended HII region associ- 
ated with Sh 252. CO line observations of Lada & Wooden 
(1979) show that this compact HII region (S252 A), the 
associated H2O maser and CO bright spot are nearly co- 
incident and located near the interface of the molecular 
cloud with the extended S252 HII region. The detection 
of the H2O maser also implies the relative youth of the 
region. Lada & Wooden (1979) also suggest S252 A to be 
in the early stages of stellar evolution in which a massive 
star (or stars) has reached the main sequence and created 
a compact HII region within its parental molecular cloud. 
VLA maps at 15 and 5 GHz of White & Gee (1986) also 
trace similar cometary morphologies. 

Using the low frequency flux densities at 1280 and 610 
MHz from our GMRT observations and 5 GHz data from 
White & Gee (1986), we derive the physical properties of 
the compact core of the HII region associated with IRAS 
06055+2039. Mezger & Henderson (1967) have shown that 
for a homogeneous and spherically symmetric core, the 
flux density can be written as 

S = 3.07 X W-^Teiy^n{l ~ e'^^""^) (3) 

where, 

riv) = 1.643a x 10^i'-^-\EM)T~^-^^ (4) 

where, S is the integrated flux density in Jy, is the elec- 
tron temperature in K, u is the frequency of observation 
in MHz, T is the optical depth, is the solid angle sub- 
tended by the source in steradians and EM is the emission 



measure in cm ^pc. a is a correction factor and we use a 
value of 0.99 (using Table 6 of Mezger & Henderson 1967) 
for the frequency range 0.6 - 5 GHz and = 8000K. The 
two GMRT maps are convolved to a common angular res- 
olution of 12" X 12" which is the resolution of the 5 GHz 
map of White & Gee (1986). In our case since the core is 
unresolved, f2 is taken as this synthesized beam size (i.e 
= 1.133 X 9a X dt, where 9a and 9), are the half power 
beam sizes). The peak flux densities of the core in the 0.6 
- 5 GHz frequency range appear to lie in the optically 
thin region. Using these peak flux densities we derive the 
emission measure for an estimated electron temperature. 

The electron temperature (Tg) of HII regions is known 
to increase linearly with Galacto-centric distance (Dg) 
(Deharveng et al 2000 and references therein). This is 
due to the decrease in heavy element abundance with Dq 
which results in higher Tg. The values of Tg derived by 
Omar et al. (2002) for a sample of three Galactic HII 
regions are also consistent with the relationship given 
in Deharveng et al. (2000). Assuming Dq as 10 kpc 
(Shirley et al. 2003) for IRAS 06055+2039 (S252A), we 
obtain a value of ~ 8000 K for the electron tempera- 
ture. Using this value of the electron temperature, the 
peak flux densities were used to fit the above equations 
(Eqns. 3 & 4). The best fit value for the emission mea- 
sure is 8.8 ± 0.4 X 10'^cm~^pc. We get an estimate of 
1.05 X lO^cm^^ for the electron density {rie — [EM/r)^-^, 
r being the core size which in this case corresponds to the 
synthesized beam size) . These values agree reasonably well 
with the estimates of Felli et al. (1977) for the brightest 
peak (A3) of the component S252 A (see Fig. 5 of FelH 
et al. 1977). They derive a value of 9.9 x 10^cm~^pc and 
5.95 X lO^cm"'^ for EM and rig respectively. They assumed 
an electron temperature of 10'*K and a distance of 2 kpc. 
It should also be noted here that the 5 GHz map of Felli 
et al. (1977) has a larger beam size (^ 8" x 21"). 

Taking the total integrated flux density of 183 mJy 
at 1280 MHz and using the formulation of Schraml & 
Mezger (1969) and the table from Panagia (1973; Table 
II), we estimate the exciting star of this HII region to be 
of the spectral type BO - BO. 5. This is consistent with the 
spectral class obtained by Felli et al. (1977) and White & 
Gee (1986). The FIR flux densities from the IRAS PSC 
yield a luminosity of ~ lO^L© which implies an exciting 
star of spectral type BO. 5, in good agreement with radio 
measurements. 

In addition to the diffuse emission seen in our 1280 
MHz map, we also detect a few discrete sources proba- 
bly representing high density clumps which are listed in 
TablelHl We designate them as Rl, R2,....R7 and their po- 
sitions are marked in Fig.|21 Three such dense clumps were 
also detected in the 5 GHz map of Felli et al. (1977). The 
position of Rl is spatially coincident with the position of 
the central bright and massive (~ BO. 5) star seen in the 
infrared cluster and is possibly the exciting source of the 
HII region. The other dense clumps could also be possi- 
ble discrete radio sources but the resolution of our map 
makes it difficult to comment on their nature. These seven 
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Table 3. Discrete sources extracted from the 1280 MHz 
map of the region associated with IRAS 06055+2039. 



Source 


RA (2000.0) 


DEC (2000.0) 


Peak Flux density 




(h m s) 


(d m s) 


(mjy/beam) 


Rl 


06 08 32.16 


+20 39 18.7 


5.7 


R2 


06 08 31.95 


+20 39 23.1 


5.4 


R3 


06 08 32.79 


+20 39 16.1 


4.3 


R4 


06 08 31.04 


+20 39 24.9 


3.3 


R5 


06 08 30.92 


+20 39 30.3 


2.8 


R6 


06 08 31.39 


+20 39 02.1 


3.2 


R7 


06 08 32.20 


+20 39 01.7 


2.8 



clumps contribute ~ 5% of the total integrated emission 
from the ionized region around IRAS 06055+2039, the re- 
maining being of diffuse nature. 

4.4. Emission from Shocked Neutral Gas 

The rotational-vibrational line of molecular hydrogen (H2 
(l-O)Sl, 2.12 /im) traces the shocked neutral gas at the in- 
terface between the ionized and the molecular gas. In the 
photo-dissociation regions (PDRs), the molecular emission 
traces the first neutral layer beyond the ionization front. 
In Fig. 1121 we compare the morphologies of the 610 MHz 
continuum map with the continuum subtracted narrow- 
band H2 (left) and Br7 (right) images. It is interesting to 
note that the shocked molecular hydrogen envelopes the 
radio emitting region. The Br7 image shows the presence 
of a faint diffuse emission which correlates well with the 
cometary morphology of the radio continuum emission. 
The H2 arc which traces the ionization front lies beyond 
the Br7 emission. Comparison with the 1280 MHz contin- 
uum map also shows similar morphology. From the posi- 
tion of the H2 arc, we estimate the radius of the HII region 
to be ~ 0.4 pc. 

4.5. Emission from Dust: Temperature, Optical Depth 
and Dust Mass 

As discussed in Sect. 14.21 the MSX images were used to 
obtain the spatial distribution of temperature and opti- 
cal depth (tio) of warm dust with the assumption that 
the dust is optically thin and the dust emissivity fol- 
lows the power law of the form eA oc A^^ (Mathis et al. 
1983; Scoville & Kwan 1976). We obtain peak values of 
1.4 X 10~^ and 155 K for tiq and the dust temperature 
respectively. Figure El shows the optical depth and the 
mid-infrared dust temperature maps. Morphologically, the 
spatial distribution of the UIBs and the optical depth con- 
tours are similar with the intensity peaks matching rather 
well. This indicates the presence of higher dust densities 
near the embedded cluster. The temperature distribution 
shows lower values near the centre with a plateau like 
feature running diagonally along the S-E and N-W direc- 
tion. Higher temperatures are towards the periphery with 
the peak seen towards the north. The optical depth and 



Table 4. Flux densities for the MSX point sources possi- 
bly associated with IRAS 06055+2039 





Flux density (Jy) 


Wavelength (fim) 


Ml' 


M2^ 


8.3 


3.52 


0.97 


12.1 


6.24 


1.40 


14.7 


2.42 


4.19 


21.3 


6.66 


13.98 



^0189. 7672+00. 3407 (0,2000 = 06''08™33.a2= ; ^2000 = 20° 39'32.04" ) 
^G189. 7677+00. 3376 (02000 = 06'' 08™ 32.40= ; ^2000 = 20°39'25.20") 



temperature are inherently anti-correlated. Hence, we see 
that the region has higher optical depth at the centre with 
the values decreasing outwards. Several peaks seen in the 
Tio map could possibly indicate the clumpy nature of the 
region. As will be evident later, we see similar trends for 
the T( 12/25) and T25 maps derived from the IRAS-HIRES 
images. 

The MSX Point Source Catalog (MSX PSC) lists two 
mid-infrared sources which fall within the radio nebulos- 
ity of IRAS 06055+2039. We designate them as Ml and 
M2. The MSX PSC flux densities for these two sources 
are listed in Table 21 Comparing the MSX mid- infrared 
colours F21/8, Fi4/i2i F14/8 and F21/14 of these two sources 
with study of the Galactic plane population by Lumsden 
et al. (2002), we infer that Ml is possibly a Herbig AeBe or 
a foreground star, whereas, M2 falls in the zone occupied 
mostly by compact HII regions. The NIR counterpart from 
2MASS catalog for M2 is the central bright IRAS point 
source. The mid- and near- infrared colours F21/8, Fg/x, 
F21/12 a-nd Fk/j of this source are also consistent with 
compact HII regions. 

The IRAS-HIRES maps (at 12, 25, 60 and 100 ^im) 
were also used to obtain the spatial distribution of warm 
and cold dust colour temperatures (T(12/25), T(60/100)) 
and optical depths {T25, tioo)- We assume the dust emissiv- 
ity to follow the power law of the form e\ oc A"-'^ . Figure [TH 
shows the dust temperature and optical depth maps. The 
maps for the optical depth T25 and colour temperature 
r(12/25) represent the warmer dust component. The dis- 
tribution is centrally dense with the optical depth peak 
and hence lower derived temperature at the centre. This 
is similar to the distribution seen in Fig. ^1 On the other 
hand rioo and T(60/100) distributions are from a rela- 
tively colder component which probably forms an envelope 
around the warmer dust. Unlike the warmer dust tem- 
perature distributions (obtained from mid-infrared emis- 
sion), the r(60/100) distribution has its peak at the centre 
with the temperature decreasing towards the periphery. 
For T25 and rioo, we obtain peak values of 2.2 x 10~^ and 
1.2 X 10^^ respectively. The dust temperature distribu- 
tions for r(12/25) and T(60/100) peak at 189 and 66 K 
respectively. Schreyer et al. (1996) derive a value of 31.2 
K for T(60/100) and 1.21 x lO^^ for noo from the IRAS 
PSC flux densities. This difference is due to the differ- 
ent resolution of the raw IRAS and the HIRES processed 
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Fig. 12. 610 MHz radio contours overlaid over the continuum subtracted H2 (left) and Br7 (right) images for the 
region associated with IRAS 06055+2039. The contour levels are same as in Fig. El The plus sign in both the images 
marks the position of the IRAS point source. 




Fig. 13. Left: The spatial distribution of dust optical depth tiq in the region around IRAS 06055+2029. The contour 
levels are at 5, 10, 20, 30, 40, 50, 60, 70, 80 and 90 % of the peak value of 1.4 x 10"". Right: The spatial distribution 
of the mid-infrared dust temperature in this region. The contour levels are for the 100,115,125,135,140,150 and 154 K. 
The peak temperature is 155 K. The plus sign in both the images marks the position of the IRAS point source. 



maps. Scaling the r25 peak value, we obtain a value of tiq 
~ 5.6 X 10~^. The difference between this value and that 
obtained from the MSX data could be a result of different 
beam sizes and/or an inhomogeneous medium. Using the 
peak value of rioo distribution, derived from the IRAS- 
HIRES maps, we estimate the warm dust mass to be 6 

Mq. 

We use the emission at submillimetre wavebands to 
study the cold dust environment in the region around 
IRAS 06055+2039. The spatial distribution of the sub- 
mm emission is shown in Fig. |3| The angular resolutions 
are 7". 8 and 15". 2 for the 450 and 850 ^m wave bands 



respectively. The main source (central dense core which 
covers the region upto ~ 25% of the peak intensity) seems 
elongated in both the maps. Apart from this dense core, 
several dust clumps are seen which could have probably 
formed due to the fragmentation of the original cloud. 

The dust mass can be estimated from the following 
relation: 

M,„,, = 1.88x10-4 (^i^^'^%^(e0 048./T._^)^2 

This is taken from Sandell (2000) and is a simpli- 
fied version of Eqn. 6 of Hildebrand (1983). The above 
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Fig. 14. Upper panel - The dust optical depth (T25) (left) and colour temperature (T(12/25)) (right) maps of the 
region around IRAS 06055+2039. The optical depth contours are 5, 10, 20, 30, 40, 50, 60, 70, 80, 90 % of the peak 
value of 2.2 x 10~^. The dust temperature contour levels are at 140, 145, 150, 155, 160, 170 and 180 K from the centre 
to the periphery. The peak value is 189 K. Lower panel - The dust optical depth (tioo) (left) and colour temperature 
(T(60/100)) (right) maps of the region around IRAS 06055+2039. The optical depth contours are 10, 20, 30, 40, 50, 
60, 70, 80, 90 % of the peak value of 1.2 x 10~^ from the centre to the periphery. The dust temperature contour levels 
are at 20, 30, 33, 39, 46, 53, 59 and 62 K. The peak value is 66 K. The plus sign in all the images marks the position 
of the IRAS point source. 



equation assumes the standard Hildebrand opacities (i.e. 
K1200GHZ = O.lcm^g^^). Here, Si, is the flux density at 
frequency i^, Td is the dust temperature which we assume 
to be 20 K (Klein et al. 2005; Mueller et al. 2002), (3 is the 
dust emissivity index and is taken to be 2 (Hildebrand 
1983) and d is the distance to the source in kpc. The 
flux densities are obtained from the JCMT-SCUBA maps 
shown in Fig. |3| To obtain the flux density of the entire 
cloud, we have integrated upto the last contour (which is 
at 5% of the peak value) . Using the above relation, we es- 
timate dust masses of ~ 70 and 90 M0 from the 450 and 
850 /im maps respectively. Assuming a gas-to-dust ratio 
of 100, the above values translate to total masses of 7000 
and 9000 Mq for the cloud from the 450 and 850 /im maps 
respectively. We also estimate the total mass of only the 



central dense core to be ~ 875 and 1250 M0 from 450 and 
850 /im maps respectively. This source (S252A) has also 
been studied by Mueller et al. (2002) at 350 /im and more 
recently by Klein et al. (2005) at 850 and 1300 /im. Scaling 
to the distance and gas-to-dust ratio assumed by us, the 
corresponding mass from Mueller et al. (2002) is ~ 600 
Mq and from Klein et al. (2005) is ^ 100 Mq. Comparing 
the mass derived from the 850 /im map, our estimate is 
very close to the mass obtained by Mueller et al. (2002) 
considering the fact that they have assumed a higher dust 
temperature (29 K). The mass estimate from the 450 /zm 
maps is 35 % lower which could have been affected by 
the large atmospheric extinction correction applied to the 
data. The mass derived by Klein et al. (2005) are rela- 
tively lower. This could be possibly due to the fact that 



14 



A. Tej et al.: Multiwavelength study of IRAS 06055+2039 




RA (JgOOO.O) 



Fig. 15. Contour maps of the emission from the ionized 
gas at 610 MHz (thick hne) and dust emission at 850 /ini 
(thin hne) are overlaid on the 2MASS Kg band image for 
the region around IRAS 06055+2039. The dashed arc rep- 
resents the position of shocked H2. The peak positions of 
the IRAS-HIRES 12 /xm (open triangle) and 100 /xm (open 
star) maps are indicated. The cross marks the position of 
the methanol and water masers. 

the flux density values presented by Klein et al. (2005) 
are lower compared to our values. From the CS line maps, 
Zinchenko et al. (1998) derive a cloud mass of 3132 Mq. It 
should be noted here that the CS maps are from a larger 
region. Furthermore, a point worth discussing here is the 
effect of varying Td and (3. Exploring the range of Td (20 
- 40 K) and /3 {1 - 2), we infer that the mass estimates 
can vary by upto a factor of ~ 8. 

4.6. Comparison of the Different Components 
Associated with IRAS 06055+2039 

In Fig. 1151 we present the various components of the re- 
gion associated with IRAS 06055+2039. The plot displays 
the contour maps of the ionized gas and the emission from 
dust overlaid on the 2MASS i^^-band image. We show 
the contour plots of 610 MHz radio emission, 850 /j,m cold 
dust emission and mark the peak positions for the 100 and 
12 iim emission from warm dust. The 850 /im emission core 
lies to the S-E of the ionized region with the warmer dust 
in between. The 12 and 100 /xm peaks lie relatively closer 
to the radio peak. The ionized region is seen to be close to 
the edge of the molecular cloud. Comparing with Fig. El 
we note that the shocked molecular gas lies in between 
the ionized region and the dense molecular core. The cen- 
tral region of the infrared cluster is located within the HII 
region which is at the edge of the molecular cloud. 

The ionized region and the dense molecular core as 
seen in the radio and the sub-mm maps respectively are 



possibly at different stages of evolution. The dense molec- 
ular core seen in our sub-mm maps does not show any 
radio emission (down to level of the rms noise which is 

0. 4 mjy/beam). The sub-mm peak is spatially offset from 
the FIR peaks and there are no MSX or NIR counterparts 
seen. This indicates a very early evolutionary stage for this 
dense and massive molecular core. It is most probably a 
very early protocluster candidate and we are sampling the 
initial collapse phase of the star forming core before the 
formation of the UCHII region (Williams et al. 2004 and 
references therein). This fact is further supported by the 
positions of the water and methanol masers which are seen 
to be coincident with the peak of the molecular core. The 
central peak of the CS map (Zinchenko et al. 1998) and 
the peaks of the 450 and 850 ^im JCMT-SCUBA maps 
match with the position of the masers. More specifically, 
it is known that methanol maser sites are generally radio 
quiet (as is the case here) and trace high mass star form- 
ing protoclusters in very early evolutionary phases (Minier 
et al. 2005). On the other hand, the ionized region which 
is associated with IRAS 06055+2039, has FIR emission, 
free-free radio emission and has NIR and MSX counter- 
parts. This could probably indicate that the source is at a 
later evolutionary stage in between an evolved protoclus- 
ter and a young cluster. Here, massive stars have started 
forming and a detectable HII region has been created. The 
cluster is partially embedded in the parental cloud. The 
sub-mm emission is weak here. A small subcluster is seen 
close to the edge of the sub-mm core spatially coincident 
with the secondary peak mentioned in Sect. I4.1TT1 Thus, 
from this multiwavelength study of the region associated 
with IRAS 06055+2039, we see the signatures of different 
evolutionary stages at different locations. 

5. Summary 

The massive star forming region associated with IRAS 
06055+2039 has been studied in detail in the infrared, 
radio and sub-mm wavelengths which lead to the following 
conclusions. 

1. The morphological details of the environment around 
IRAS 06055+2039 show that we have probed different 
stages of evolution of star formation present in this 
cloud complex. 

2. High sensitivity and high resolution radio continuum 
maps at 1280 and 610 MHz obtained from our observa- 
tions using GMRT show interesting cometary morphol- 
ogy. Apart from the diffuse emission, the 1280 MHz 
continuum map also shows the presence of several dis- 
crete sources which probably represent high density 
clumps. The total integrated emission implies an ex- 
citing star of spectral type BO - BO. 5 which is consis- 
tent with the estimates derived from NIR colours and 
IRAS PSC flux densities. 

3 The dense molecular cold core has been probed with 
JCMT-SCUBA at 450 and 850 /im. The sub-mm peak 
is spatially offset from the peak of radio emission. The 
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sub-mm emission is most probably from a very early 
protocluster candidate and the radio emission is from a 
region at a later evolutionary stage where the massive 
stars have formed a HII region though the cluster is 
still partially embedded in the prenatal cloud. 

4. In the NIR, the Br7 emission correlates well with 
the radio continuum emission. The H2 (1-0) SI line 
of molecular hydrogen, which traces the first shocked 
neutral layer beyond the ionization front is scon as an 
arc towards the N-E of the central IRAS point source 
and envelopes the ionized emission mapped at radio 
and NIR wavelengths. The shocked molecular gas lies 
between the ionized region and the dense cloud core. 

5. Using the 2MASS data, we derive a power law slope of 
0.43±0.09 for the KLF of the NIR cluster associated 
with this star forming complex. This is consistent with 
the values obtained for other young embedded clusters. 
Wc estimate an age of 2 - 3 Myr for the cluster. The 
physical structure of the cluster suggests that it is not 
yet in complete dynamical equilibrium which is consis- 
tent with the scenario of different evolutionary stages 
seen in the complex. 

6. The spatial distribution of the emission from the 
UIBs as extracted from the MSX images also displays 
cometary morphology. The total UIB emission from 
MSX band A is estimated to be 2.85 x IQ-'^'^Wm.-'^ 
and is ^ 39% of the emission estimated from the IRAS 
LRS spectrum. The MSX mid-infrared colours of the 
central bright IRAS point source corresponds to a com- 
pact HII region. 

7. The spatial distribution of temperature and optical 
depth of interstellar dust has been presented based on 
the mid- and far- infrared data from the MSX and 
IRAS (HIRES) surveys. From the derived peak value 
of the Tioo distribution, we estimate the warm dust 
mass to be ^ 6 . From the sub-mm emission at 450 
and 850 /xm, we estimate the total mass of cloud to be 
~ 7000 - 9000 Mq. 
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